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We use numerical simulations of different dark energy cosmologies to investigate the 
concentration-mass (c — M) relation in galaxy clusters. In particular, we consider a reference 
A cold dark matter (ACDM) model, two models with dynamical dark energy, viewed as a 
quintessence scalar field [using a Ratra and Peebles (RP) and a supergravity (SUGRA) potential 
form], and two extended quintessence models, one with positive and one with negative coupling 
(EQp and EQn respectively), where the quintessence scalar field interacts non-minimally with 
gravity (scalar-tensor theories). All the models are normalized in order to match CMB data from 
Wilkinson Microwave Anisotropy Probe 3 (WMAP3). For each model, we have performed nu- 
merical simulations in a cosmological box of (300 Mpc /i -1 ) 3 . We fit the dark matter profile 
with a Navarro-Frenk- White (NFW) profile, and recover the concentration of each halo. We con- 
sider both the complete catalog of clusters and groups and a subsample of relaxed objects. The 
c — M relation of our reference ACDM model is in good agreement with the results in literature, 
and relaxed objects have a higher normalization and a shallower slope with respect to the com- 
plete sample. For the different dark energy models, we find that for ACDM, RP and SUGRA 
the normalization of the c—M relation is linked to the growth factor D + and the power spectrum 
normalization Og, with models having a higher value of <J%D + having also a higher normalization. 
This simple scheme is no longer valid for EQp and EQn because in these models it is present a 
time dependent effective gravitational interaction, whose redshift evolution depends on the sign 
of the coupling. This leads to a decrease (increase) of the expected normalization in the EQp 
(EQn) model. This result shows a direct manifestation of the coupling between gravity and the 
quintessence scalar field characterizing EQ models that cannot be seen at the background level 
but can be investigated in the non-linear regime. 
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1. Introduction 

Over the last decade great observational evidence ([jl]-01) has shown that at the present time the 
Universe is expanding at an accelerated rate. This fact can be attributed to a component with neg- 
ative pressure, which is usually referred to as dark energy, that today accounts for about 3/4 of the 
entire energy budget of the Universe. The simplest form of dark energy is a cosmological con- 
stant term A in Einstein's equation, within the so-called A cold dark matter (ACDM) cosmologies. 
Though in good agreement with observations, a cosmological constant is theoretically difficult to 
understand in view of the fine-tuning and coincidence problems. A valid alternative consists in a 
dynamical dark energy contribution that changes in time and space, often associated to a scalar field 
('quintessence') evolving in a suitable potential ([Q, |6|]). Dynamical dark energy allows for appeal- 
ing scenarios in which the scalar field is the mediator of a fifth force, either within scalar-tensor 
theories or in interacting scenarios ([0-|lJ and references therein). In view of future observations, 
it is of fundamental interest to investigate whether dark energy leaves some imprints in structure 
formation, giving a practical way to distinguish among different cosmologies ([|2|-[F7|]). 

The internal properties of dark matter halos are known to reflect their formation history and 



thus the evolution of the background cosmology. Q18| ] (hereafter NFW) found that the dark matter 
profile of a halo can be characterized by a scale radius, which is linked to the virial radius through 
the concentration of the object. The concentration of a dark matter halo is related to the mean 
density of the universe at the halo formation time. 

Because of the hierarchical nature of structure formation and the fact that collapsed objects 
retain information on the background average matter density at the time of their formation ([p"8|]), 
concentration and mass of a dark matter halo are related. Since low-mass objects form earlier 
than high-mass ones, and since in the past the background average matter density was higher, 
low-mass halos are expected to have a higher concentration compared to high-mass ones. These 
expectations have been confirmed by the results of Af-body numerical simulations which find, at 
z = 0, a concentration-mass relation c(M) °c M a , with a ~ —0.1 ([|T^, [20| , 0]), with a log-normal 



scatter ranging from 0.15 for relaxed systems to 0.30 for disturbed ones ([21 ]). 

A^-body simulations have been carried out by several authors in order to study the c — M re- 
lation in dark matter halos with sizes of galaxy groups and clusters. [|^] performed simulations 
with different cosmological models in order to verify the effects of dark energy dynamics. They 
found that the halo concentration depends on the dark energy equation of state through the linear 



growth factor at the cluster formation redshift, D + {z co u). Ij22p also noted that non-relaxed objects 
have a lower concentration and a higher scatter with respect to relaxed ones. [ |23| ] made a compari- 
son between concentrations in the Wilkinson Microwave Anisotropy Probe 1 (WMAP1), WMAP3 
and WMAP5 cosmologies in order to study the effects of different cosmological parameters (in 
particular the power spectrum normalization as) on the c — M relation. 

Since the concentration of a halo is linked to the background density of the universe at the 
time it collapsed, and since different dark energy models predict different evolutions of the cosmo- 
logical background, it is interesting to investigate the impact of dark energy on the c — M relation. 
Moreover, since some dark energy models can also affect the linear and non linear evolution of the 
density fluctuations, leaving some imprints in collapsed structure, one can think about using the 
c — M relation as a cosmological probe, orthogonal to others that are commonly used, jty-body cos- 
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mological simulations of extended quintessence models, including the effects on the c — M relation, 
were presented in [|4j] . 

2. Dark energy models 



We consider the same cosmological models discussed in [25]. Here we recall only the main features 



of the different models, and refer to [25 1 for more details. 

As a reference model we use the concordance ACDM model. This model is characterized by 
the presence of a dark energy component given by a cosmological constant A, with equation of 
state wa = — 1 ■ 

The second case is a model with dynamical dark energy, given by a quintessence scalar field 
(j) with an equation of state evolving with redshift, w = w(z) ([§, ||]). As in [25], as potentials for 



minimally coupled quintessence models, we consider an inverse power-law potential 

M 4+a 

V(0) = — , (2.1) 
the so called RP potential ([0]), as well as its generalization suggested by supergravity arguments 



([26]), known as SUGRA potential, given by 

M A+a 

V($) = -^ r exp(4nG<j) 2 ) , (2.2) 

where in both cases M and a > are free parameters (see Table [T]for details). 

The third possibility we consider is the case in which </> interacts non minimally with gravity 
'lH- Ell)- m particular we refer to the extended quintessence (EQ) models described in [27], [28] 



and [jlO|]. The parameter ^ represents the "strength" of the coupling (see Table |]for details). In 
particular we consider here a model with positive coupling <^ > (EQp) and one with negative 
£, < (EQn). For an extensive linear treatment of EQ models we refer to [|10|]. Here we only 
recall for convenience that EQ models behave like minimally coupled quintessence theories in 
which, however, a time dependent effective gravitational interaction is present. In particular, in the 
Newtonian limit, the gravitational parameter is redefined as 

2[F + 2(dF/d^) 2 } 1 

[2F + 3(dF/d$) 2 } SkF ' ^ ; 

Here the coupling F((f>) is chosen to be 

F(0) = i + <^ 2 -«/> 2 ), (2.4) 



with K = 87TG*, where G* represents the "bare" gravitational constant ([[29|]). 



For small values of the coupling, that is to say <C 1, the latter expression becomes 

-^-~l-87TG^(4> 2 -4> 2 ) , (2.5) 

which manifestly depends on the sign of the coupling ^ . We note that, since the derivative of the 
RP potential in equation (2.1) with respect to (f> is dV{ty)/dty < 0, we have (j) 2 < 0q. This leads to 



the behaviour of G/G* shown in Fig. [j]. Note that the corrections are only within the percent level. 
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Figure 1: Correction to the gravity constant for the two extended quintessence models, EQp (cyan) and EQn 
(red), as expressed in equation (2.5). Note that the corrections are only within the percent level. 

3. Numerical simulations 

In order to study the formation and evolution of large scale structures in these different cosmologi- 
cal scenarios we use Af-body simulations performed with the GADGET-3 code ([^0|, |]]]). For each 
model, we simulated a cosmological box of size (300 Mpc h~ 1 ) 3 , resolved with (768) 3 dark matter 
particles with a mass of m& m « 4.4 x 10 9 M Q h^ 1 . 

As in [p^], we modified the initial conditions for the different dark energy scenarios adapting 
the initial redshift for the initial conditions in the dark energy scenarios determined by the ratio of 
the linear growth factors D + (z), 

D+{Zini) _ £>+,ACDm(Zacdm) ,3 ^ 

D+(0) D +iACD M(0) 
Therefore, all simulations start from the same random phases, but the amplitude of the initial fluc- 
tuations is rescaled to satisfy the constraints given by CMB. 

Our reference ACDM model is adapted to the WMAP3 values ([{3^]), with the following cosmo- 
logical parameters: 

• matter density: £2om = 0.268 

• dark energy density: ^oa = 0.732 

• baryon density: Cloy = 0.044 

• Hubble parameter: h = 0.704 

• power spectrum normalization: 0% = 0.776 

• spectral index: n s = 0.947 

We trimmed the parameters of the four dynamical dark energy models so that wq = w(0) « 
—0.9 is the highest value still consistent with observational constraints in order to amplify the 
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Model 


a 




wo 




ACDM 






-1.0 


0.776 


RP 


0.347 




-0.9 


0.746 


SUGRA 


2.259 




-0.9 


0.686 


EQp 


0.229 


+0.085 


-0.9 


0.748 


EQn 


0.435 


-0.072 


-0.9 


0.729 



Table 1: Parameters for the different cosmological models: a is the exponent of the inverse power-law 
potential; £ is the coupling in the extended quintessence models; wo is the present value of the equation of 
state parameter for dark energy; <J% is the normalization of the power spectrum. 
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Figure 2: Redshift evolution of the equation of state parameter w for the different cosmological models 
considered: ACDM (black), RP (blue), SUGRA (green), EQp (cyan), and EQn (red). 

effects of dark energy. Fig. || shows the evolution with redshift of w in each cosmology. The 
parameters £2om, &oa, &0b, h, and n s are the same for all the models, but since we normalize all 
the dark energy models to CMB data from WMAP3, this leads to different values of as for the 
different cosmologies: 

<D+,acdm(zcmb) ,~ ~s 

<*8,DE — 08.ACDM— r; j r— , (P-A) 

^+,de(Zcmb) 

assuming zcmb = 1089. This fact, along with the different evolution of the growth factor D + 
(shown in Fig. |5|), has an impact on structure formation. Table [l] lists the parameters chosen for the 
different cosmological models. 

Using the outputs of the simulations, we extract galaxy clusters from the cosmological boxes, 
using the spherical overdensity criterion to define the collapsed structures. We take as halo centre 
the position of the most bound particle. Around this particle, we construct spherical shells of 
matter and stop when the overdensity drops below 200 times the mean (as opposed to critical) 
background density defined by Q m po c ; the radius so defined is denoted with /?200m and the mass 
enclosed in it as Mzoom- We consider all the halos having M200m > 10 14 M Q h . In addition, we 
selected subsamples of the 200 objects with M2oom closest to 7 x 10 13 M Q h~ v , 5 x 10 13 M hr l , 
3 x 10 13 M h , and 10 13 M Q h . Starting from the centres of the halos, we construct radial 
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Figure 3: Ratio between the value of <J%D + for the ACDM (black), RP (blue), EQp (cyan), EQn (red), and 
SUGRA (green) cosmologies and the corresponding value for ACDM as a function of redshift. 



profiles by binning the particles in radial bins. We concentrate on objects at z = 0. For the following 
analysis, we also calculate for each cluster selected in this way the radius at which the overdensity 
drops below 200 times the critical background density and denote it as /?200- The corresponding 
mass is indicated as M200. It is useful to define a quantitative criterion to decide whether a cluster 
can be considered relaxed or not because, in general, relaxed clusters have more spherical shapes, 
better defined centres and thus are more representative of the self-similar behaviour of the dark 
matter halos. We use a simple criterion similar to the one introduced in [22]: first we define x // as 
the distance between the centre of the halo (given by the most bound particle) and the barycentre of 
the region included in /?200 m ; then we define as relaxed the halos for which x ff < 0.07/?200m- We 
plot the distribution of x (/ ff for the objects in the five cosmological models at z = in the left-hand 
panel of Fig. |] Note that the distribution and the median value of x (/ ff are similar in the different 
cosmological models. 



4. c — M relation 

For each cluster at z = in the five cosmological models under investigation, we perform a log- 
arithmic fit, using Poissonian errors (In 10 x ^/n^)^ 1 (where n^m is the number of dark matter 
particles in each radial bin, of the order of 10 — 10 3 depending on the mass of the object), of the 
three-dimensional dark matter profile p c i m (r) in the region [0.1 — l]/?200 (where the value of /?2oo is 



taken directly from the true mass profile) with a NFW profile ([|18|]) 

Pdmjr) 8 

p c ~(r/r,)(l + r/r,) 2 ' 



(4.1) 



where p c is the critical density, r s is the scale radius and 8 is a characteristic density contrast. 
Then, instead of defining C200 = ^200 A s . we directly find the concentration parameter C200 from the 
normalization of the NFW profile 
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Figure 4: Left-hand panel: the distribution of x /f (in units of ^200m) f° r the objects in ACDM (black), RP 
(blue), SUGRA (green), EQp (cyan), and EQn (red) at z = 0. The vertical lines of the corresponding colours 
mark the median value of x a ff in each cosmological model. The vertical pink line corresponds to the value 
defining relaxed objects, x ff = 0.07/?200m- Right-hand panel: the same as in the left-hand panel, but for 



200 c 



3 

200 



C200 



(4.2) 



[ln(l + c 20 o) l+C2m _ 

We require the central density parameter 8 to be greater than 100 and the scale radius r s to be 
within [0.1 — l]/?200- We exclude the inner regions from the fit because we are limited in resolution 
inside a given radius. We indicate the dark matter concentration found in this way as C20Odm- We 
define the rms deviation o rms as 

J Nbins 

°lm S = Tt E [toglOPdm; ~ logloPiVFW,] 2 , (4-3) 

Nbins i=\ 

where Nbins is the number of radial bins over which the fit is performed and Pnfw is the best-fitting 
NFW profile. We plot the distribution of o rms for the objects in the five cosmological models at 
z = in the right-ended panel of Fig. |] Note that the distribution and the median value of o rms are 
similar in the different cosmological models, meaning that the NFW profile is as good as in ACDM 
in describing the dark matter profile of galaxy clusters in dark energy cosmologies. 

We bin the objects in the complete sample in groups of 200, so that we have bins around 
10 13 M© r',3x 10 13 M Q r',5x 10 13 M Q hr l , and 7 x 10 13 M h~ l . For halos more massive 
than 10 14 M h~ l , we bin the objects starting from the low-mass ones, so that the most massive 
bin can contain less than 200 objects. The analysis for the relaxed sample is done by selecting the 
relaxed objects inside each bin. Once we have c^mdm f° r eac h object in each mass bin, since the 
distribution of C2oodm is log-normal inside each bin, we evaluate the mean M200 an d the mean and 
rms deviation of logioC2oo<f»i in each bin, for the two samples. In the following of the proceeding, 
when we indicate the value of C2oodm in a mass bin, we refer to iO< lo gi° c 2ooA»>. 

With the mean and rms deviation of logioC2oorfm in each bin at hand, we fit, for the complete 
and relaxed samples, the binned c — M relation using 
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logio^oo = logio^ + B logio ^ Ypn^ ) ' < - 4 " 4) 

where logioC2oo and M200 are the mean values in each bin. For the error on the mean of log 10C200A?! 
in each bin, o E , we use the rms deviation of logioC2oo^m divided by the square root of the number 
of objects in the bin. For each fit we also define 

N^as S /logi C200rfm, - logl0C200/if, \ 2 . _ 

X = h ( % ) • W 51 

where N mass is the number of mass bins over which the fit is performed and C2oofit is obtained 
from the best fit of equation (|4.4[), and evaluate the reduced chi-squared % 2 , i.e. % 2 divided by the 



number of degrees of freedom. 

In the reference ACDM model, relaxed objects have a higher normalization and a shallower 
slope with respect to the complete sample (see Table || for details). By comparing our results for 



ACDM with previous works in literature (see [ |33| ] for the comparison), we find a good agreement, 



in particular when the values of the cosmological parameters are similar, as in [p3[]. Thus, when 
comparing the impact of different dark energy models on the c — M relation, we can rely on our 
ACDM model as a reference. The c — M relation for galaxy clusters extracted from dark matter 
only simulations of different dark energy models, including RP and SUGRA, has been studied in 
Jl9|]. They fit a formula similar to equation ([T4]) and find that, when the same as is used for all the 
models, the normalization of the c — M relation for dark energy cosmologies is higher compared to 
ACDM, depending on the ratio between the growth factors through 

A ^ A D +,De{ZcoIi) „ 

^DE r ^Iacdm t: -, r, (4.0) 

L>+,ACDM\Zcoll) 



where the collapse redshifts z co ll are evaluated following the prescriptions of [J34|]. When as values 
are normalized to CMB data, as we do in this work, the normalization of the c — M relation for dark 
energy cosmologies is lower compared to ACDM. We find that, in order to recover the values of 



the normalization they quote in this case, equation ( |4.6| ) should be multiplied by the ratio between 
the values of as, i.e. as.DE/^s.ACDM- This fact goes in the same direction as what found in [|23|], 
where models with higher as also have a higher normalization of the c — M relation. 



5. Results 



We compare the c — M relation for the dark energy models under investigation with the one derived 
for the ACDM cosmology. In Table || we summarize the best-fitting parameters, the standard 
errors and the reduced chi-squared of the c — M relation equation ([k4]) for the five cosmological 
models here considered, both for the complete and relaxed samples. For the complete sample, the 
differences in the normalization A between ACDM and the other cosmological models are less than 
10%, with EQn being the only model having a higher normalization. The slope B is within 5% of 
the ACDM value for all the models with the exception of EQn, which shows a 30% flatter slope. For 
the slope the differences among the models, excluding EQn, are smaller than the standard errors, 
while for the normalization these differences are significant. If we limit ourselves to the best-fitting 
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Model ff 8 


A a A 


B a B 


X 2 


A a A 


B o B 


X 2 




all 


relaxed 


ACDM 0.776 
RP 0.746 

SUGRA 0.686 
EQp 0.748 
EQn 0.726 


3.59 0.05 
3.54 0.05 
3.41 0.05 
3.36 0.05 
3.70 0.05 


-0.099 0.011 
-0.103 0.011 
-0.098 0.013 
-0.097 0.012 
-0.069 0.013 


0.48 
1.14 
1.50 
0.35 
0.78 


4.09 0.05 
4.08 0.05 
3.94 0.06 
3.84 0.05 
4.25 0.06 


-0.092 0.011 
-0.081 0.011 
-0.081 0.012 
-0.097 0.011 
-0.081 0.013 


0.66 
0.92 
1.55 
1.32 
0.51 



Table 2: Best-fitting parameters, standard errors and reduced chi-squared % 2 of the c — M relation equation 
(4.4) for dark matter density profile fit in the region [0.1 — 1]^200 for the complete and relaxed samples of 
the five different cosmological models at z — 0. 



values, given that the slope is almost identical and that all the cosmological parameters except o% 
are fixed, we expect that the normalization should follow the values of Og, i.e. the higher 0$ the 
higher the normalization (see [|3|]), and D + , i.e. the higher D + at z co ll the higher the normalization 
(see Jl9|]). The quantity controlling the normalization is thus expected to be o%D + (z C ou), which is 
plotted as a function of redshift in Fig. || for the five cosmological models. Independently of the 
precise definition of z co iu the cosmological model with the highest value of this quantity is ACDM, 
followed by RP, EQp, EQn, and SUGRA. We do expect the normalization of the c — M relation to 
follow the same order, with ACDM having the highest and SUGRA the lowest. Instead we see that, 
on the one hand, EQp which has the third highest o%D + has the lowest normalization while, on the 
other hand, EQn which has the second lowest C%D + has the highest normalization. The relative 
order of o%D + and A is preserved for ACDM, RP and SUGRA, as in JEfl. 

For the relaxed sample, compared to ACDM, the differences in the normalization are less than 10%, 
while the differences in the slope can almost reach 15%, but they are compatible with the standard 
errors. Also in this case, the most extreme cosmologies are EQp and EQn, whose normalization 
goes in the opposite direction with respect to their C%D + . This fact confirms the conclusions we 
have drawn from the complete sample. The values of the reduced chi-squared indicate that equation 
(4.4) is a good parametrization of the c — M relation for almost all cosmological models. Only 
SUGRA has high values both for the complete and relaxed samples. 

Our results are in good qualitative agreement with the findings of [24], where halos in extended 
quintessence models have lower (higher) concentrations with respect to the ACDM case for positive 
(negative) values of the scalar field coupling. 

We plot the best-fitting c — M relation for all the cosmological models, along with the binned 
data, in Fig. |5[ We clearly see that the results on the normalization are due to differences in the 
concentrations over a wide mass range. If we look, for example, at the complete sample (left-hand 
panel of Fig. |J), we see that the different slope of EQn is mainly originated by the less massive 
bin. But with the exception of this bin, EQn shows the highest concentration in almost all the mass 
bins, while in general EQp has the lowest concentration. For the relaxed sample (right-hand panel 
of Fig. |5|), the relative behaviour of the different cosmological models is even clearer, and indeed 
the differences in the slope are less pronounced. 

Before drawing our conclusions about the EQ models, we want to take into account the de- 
pendence of the normalization on the slope that characterizes the c — M relation in the different 
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all z = relaxed z = 




10 13 10 u u 10' 3 10'" 

M 200 ( M G>) M 200 ( M o) 



Figure 5: Left-hand panel: the values of C2oodm f° r the complete sample of the ACDM (black), RP (blue), 
SUGRA (green), EQp (cyan), and EQn (red) cosmologies at z = 0. The lines of the corresponding colours 
are our best fit of c — M relation equation (4.4) and the vertical black bar is the error on the normalization 
of ACDM as listed in Table 2. The symbols in the low part of the panel are the ratios between cjoodm f° r 
the model and C2oodm f° r ACDM. Right-hand panel: the same as in the left-hand panel, but for the relaxed 
sample. 



cosmological models. To do this, we fix the slope at the best-fitting value for the complete sample 
of ACDM at z = (i.e. B = -0.099, see Table g) and we fit equation ( p~4| ) with only A as a free 
parameter. We report the results in Table [3] and plot them in Fig. [6], which summarizes almost all 
the information on the c — M relation at z = for the cosmological models under investigation. We 
show the values of the reduced chi-squared of the fit as a reference, but we do not discuss them 
because we are imposing the slope for ACDM also to other models. Also in this case, relaxed 
objects have a higher normalization compared to the complete sample. Then, as a general trend, 
both fixing or keeping the slope free, the normalization is decreasing going from ACDM to RP to 
SUGRA, independently of the dynamical state. Finally EQn always has the highest normalization 
while EQp alway has the lowest. The behaviour of ACDM, RP and SUGRA is in agreement with 
the simple idea that the normalization of the c — M relation is driven by the value of o&D + , but the 
one of EQp and EQn is not. 

We hint that the behaviour of EQp and EQn is linked to the redshift evolution of the effective 
gravitational interaction G, as pointed out in Section ||. In fact, in contrast with ACDM, RP and 
SUGRA, in EQ models the gravitational constant G is substituted by G, which is higher (lower) 
than G at high redshift for positive (negative) values of the coupling constant B, , while it is equal 
to G at z = in order to recover General Relativity (see Fig. [I]). This means that in EQp gravity 
becomes weaker at low redshift compared to high redshift, while in EQn it becomes stronger. Thus 
one can expect that in EQp (EQn) the halos expand (contract) due to the change in the gravitational 
interaction, resulting in lower (higher) concentrations with respect to the case in which gravity is 
constant. 



6. Conclusions 

In this proceeding, we reviewed the c — M relation for the halos extracted from the simulation set 
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Model 


OS 


A 


Ok 


B 


o B 


X 2 


A 


Ok 


B o B 


X 2 




all 


relaxed 


ACDM 


0.776 


3.59 


0.05 


-0.99 


0.11 


0.48 


4.08 


0.04 


-0.99 — 


0.62 


RP 


0.746 


3.55 


0.04 


-0.99 




0.97 


4.04 


0.05 


-0.99 — 


1.19 


SUGRA 


0.686 


3.41 


0.04 


-0.99 




1.20 


3.89 


0.04 


-0.99 — 


1.69 


EQp 


0.748 


3.36 


0.04 


-0.99 




0.30 


3.83 


0.04 


-0.99 — 


1.11 


EQn 


0.726 


3.65 


0.04 


-0.99 




1.56 


4.21 


0.05 


-0.99 — 


0.74 



Table 3: Best-fitting parameters, standard errors and reduced chi-squared % 2 of the c — M relation equation 
(4.4), with B fixed at the best-fitting value for the complete sample of ACDM at z = 0, for dark matter density 
profile fit in the region [0.1 — l]/?200 for the complete and relaxed samples of the five different cosmological 
models at z = 0. 
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Figure 6: Left-hand panel: best-fitting normalization comparison for equation (4.4) for the ACDM (black), 
RP (blue), SUGRA (green), EQp (cyan), and EQn (red) cosmologies. Triangles: dark matter profile fit, 
complete sample. Squares: dark matter profile fit, relaxed sample. The vertical black bar is the error on the 
normalization of the complete sample of ACDM. Right-hand panel: the same as left-hand panel but with B 
fixed at the best-fitting value for the complete sample of ACDM at z = 0. 



introduced in [25] and [33]. We find that the normalization of the c — M relation in dynamical dark 
energy cosmologies is different with respect to the ACDM one, while the slope is more compatible. 
In particular, at z = 0, the differences in the normalization for RP and SUGRA when compared to 
ACDM reflect the differences in C%D + , with models having a higher OgD + also having a higher 
normalization. This simple scheme is not valid for the EQp and EQn scenarios. In the former case, 
the normalization is lower than expected considering C%D + , while in the latter it is higher, and 
indeed EQn is always the model with the highest normalization, regardless of the dynamical state 
of the objects. This behaviour is due to the different redshift evolution of the effective gravitational 
interaction G that characterizes these models. Indeed, going from high to low redshift, G decreases 
(increases) for EQp (EQn), making the halos expanding (contracting) and thus decreasing (increas- 
ing) the concentration. This is a very important result because it shows a direct manifestation of the 
coupling between gravity and the quintessence scalar field that cannot be seen at the background 
level but can be investigated in the non-linear regime. 
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